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We studied the semiconductor response with respect to high intensity resonant excitation on 
short time scale when the contribution of the Fermi statistics of the electrons and holes prevails. 
We studied both the single and double pulse excitations. For the latter case we considered the time 
evolution of the multi-wave mixing exciton polarization. The main difference between the excitation 
by a single pulse or by two non-collinear pulses is that the Rabi oscillations of the multi-wave mixing 
response are characterized by two harmonics. Analyzing the operator dynamics governed by the 
external excitation we found that there are three invariant spin classes, which do not mix with the 
evolution of the system. Two classes correspond to the bright exciton states and one contains all 
dark states. We found that the dynamics of the classes is described by six frequencies and the 
Rabi frequencies are only two of them (one for each bright class). We discuss the effect of the 
dispersion of the electrons and holes and the Coulomb interaction describing the semiconductor by 
the semiconductor Bloch equation (SBE). We show that if initially the system is in the ground state 
then the SBE preserves the invariant spin classes thus proving absence of the dark excitons in the 
framework of this description. We found that due to the mass difference between holes of different 
kind additional Rabi frequencies, two of those present in the operator dynamics, should appear in 
the evolution of the exciton polarization. 

PACS numbers: 71.35.-y,71.45.Gm,78.47.Fg 



One of the main tools of probing the complex character of the many-body correlations and interactions in semicon- 
ductors is the multi-wave mixing response. The multi-wave mixing polarization in optically excited semiconductors 
produces the signal in directions that are prohibited in the linear regime because of momentum conservation, thus giv- 
ing access to the semiconductor many-body excitations. The typical example is the four-wave mixing spectroscopy^^ 
Recently also wave mixing of higher orders started to draw attentioni^i^i^ A great success in understanding the 
mechanism of the formation of the nonlinear response in general was achieved in relatively low-field limit using the 
perturbation theory with respect to the external field. -"ii^i^iiSiiiii^ Perturbational description, however, is not suitable 
for investigating the coherent reconstruction of the spectrum, such as in the case of Rabi oscillations. For a qualita- 
tive analysis of such phenomena few-level quantum models have been usedi ^^'^'^i^^ This approach, however, misses the 
important property of the excitations in semiconductors. These excitations constitute a quantum field rather than a 
canonical quantum mechanical system. The system supports infinitely many states unless, of course, the spectrum 
has truly discrete component, that is when the excitations are localized, such as in quantum dots. Thus a description 
non-perturbative with respect to the external field must deal with many-body aspects of the dynamics of the systemi^ 
In order to treat this problem a variety of methods based on derivation of respective closed equations of motion was 
de^^e\o•ped^ ^ ' ^^^ ' ^'^ ^^^- ' ^^■ ' ^^^'^^ Unfortunately, the dynamical equations turn out to be very complex owing to the Coulomb 
interaction, which is shown to be crucially important at relatively long time scales<22i2^ Therefore, one has to resort 
to numerical calculations investigating the dynamics of the exciton polarization. The numerical simulations proved 
to describe successfully the dynamic o^'^'^^i^'* but they are difficult to apply for studying the detailed effect of different 
contributions to the dynamics and its dependence on parameters of the system. As a result, the theory of the semi- 
conductor response suffers from a lack of exact results obtained in controllable approximations which would guide the 
respective numerical, theoretical, and experimental studies. 

At short time scales, however, one can rely on significant simplification of the dynamics of the system due to 
negligible phase change during the optical excitation. Effectively, the system follows the external field, which is 
illustrated best by a simple dynamical model 



with uj{P) being a real positive function. If initially the system is at rest, P(0) = 0, we can present the time 
dependence P{t) in the form 
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This representation is convenient for comparing the responses of the system with respect to the ^-functional excitation 
E{t) — eSit + 0) and to the excitation with the piece-wise constant ampUtude E{t) — e/T, where T is the duration of 
the forced regime. As follows from Eq. ([2]), in the case of the 5-pulse immediately after the excitation is switched off 
one has P{t) — e. For the excitation with constant amplitude Eq. Q yields an estimate maximally different from the 
(5-pulse case in the form P{T) = e(l - e'^'^^^) jiijj^xT , where LOm = supo<p<g ti^(P). Thus if the excitation pulse is 
shorter than the time scale determined by Wm, i.e. if LOmT ^ 1, then the response with respect to the constant pulse 
of finite duration differs up to quadratic terms from the (5-pulse response only by the phase factor e^'"'"^/^. 

An immediate application of this consideration to the evolution of the exciton polarization in optically excited 
semiconductors is prevented by the more complex character of the polarization dynamics, namely, by the field induced 
coupling between the polarization and the charge densities. In order to see the principal differences introduced by 
this coupling we consider a model 

P(t) = -liuPit) - iE{t)n{t) - iE{t), 

h{t) = mn{t) - iE{t)P{t). ^ ' 

Straightforward substitution of the ^-shaped pulse into this equation obviously fails since now the amplitudes in front 
of the (5-function in the coupling terms are given by the quantities, n{t) and P{t), which experience discontinuity at 
exactly the point of singularity of the (5-function. Therefore, we need to consider the case of the external pulse of 
finite duration. For the pulse with piece- wise constant amplitude E = e/T, as in the previous example, Eq. ([3]) with 
the initial conditions P(0) = 0, n(0) = is solved for t <T hy 
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P(t) = - i— sin(m/2) sm{Qt/2) - i- cos(m/2) 



(4) 



where ft^ = uj"^ /A+E"^. It is seen that in the limit of the excitation pulses short compared to the typical dynamical time 
scale, ujT <C 1, one has to distinguish the cases of weak and strong excitations, ojT/e > 1 and ivT/e <C 1 respectively. 
In the first case up to the quadratic terms, fi^T^, the response P{T) looks like the one of the system where P and 
n are decoupled with respect to (5-pulse excitation. This illustrates the relation with the 5-functional approximation 
popular in studying the nonlinear optical response by using perturbational approach, which corresponds to weak 
scattering. 

In the limit of strong external field when ujT ^ 1 while e > I the evolution of the system is significantly different. 
Neglecting the terms quadratic in LuT/e we obtain that the evolution of the system is similar to the one we have in 
the case lo — (i.e. Rabi flopping) with the only difference that P{t) acquires a (small) imaginary part. Applying 
now the similar idea for estimating a phase shift as for the analysis of Eq. ([T]) we can see that the conclusion about 
small perturbations vanishing with T — > holds also in the case of P-dependent frequency and renormalized coupling 
E ^ E + r]{\P\,\n\). 

These results serve as a general background for the analysis of the basic features of the immediate response of a 
semiconductor excited by short pulses of high intensity. In this regime, while taking exactly into account the many- 
body effects due to the fermion nature of the electrons and holes, one can neglect the Coulomb interaction and still 
be able to provide a qualitative description. From this perspective the roles played by statistics and interactions are 
clearly different. The statistics impose instantaneous constraints on the dynamics. Interactions, in turn, need time to 
develop their effect. 

Explicitly the possibility to neglect the Coulomb interaction in the semiconductor Bloch equation (SBE) was shown 
for the two-band model in Ref. 25. In Ref. ,2^ the numerical investigation of the coupled Maxwell-Bloch equations 
in multiple quantum well structures has shown that one can reproduce with good accuracy a few Rabi flops at short 
time scale. The case of very short excitation of very high intensity was studied in Ref. [2^ from the perspective 
of the carrier-wave Rabi flopping ^Ii^iSi It was demonstrated that even the spectrum of the emitted radiation is 
approximated remarkably well by the free-carrier model, which corresponds to neglecting the Coulomb interaction in 
the SBE. It should be noted, however, that this extreme limit should be considered with care because of the failure 
of the rotating wave approximation. ?-'3'^ 

In the present paper we extend the consideration of the semiconductor response with respect to short intensive 
excitation taking into account multiple hole states in the valence band and the formation of the multi-wave mixing 
response. The structure of the paper is as follows. In Section [II] we consider the limiting case when the effect of the 
internal dynamics can be completely neglected, which corresponds to the cases $7 = in the examples considered 
above. We obtain classes of state related to the spin selection rules and related to these classes the spectrum of the 
Rabi frequencies. In Section [IIII we discuss the effect of the internal dynamics and the Coulomb interaction. 
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II. RABI OSCILLATIONS 



We consider the dynamics of excitation of a semiconductor quantum well during the action of an ultrashort pulse 
of high intensity so that the period of the Rabi oscillations is smaller than the typical time scales determined by 
the internal dynamics. We consider the question of validity of this approximation in Section IIIII We assume that 
the pulse is tuned to resonance with the exciton levels lying below the band edge and that the envelope of the 
pulse is constant while the field is on. By means of a canonical transformation the harmonic time dependence of 
the excitation field can be excluded, (see e.g. Ref. Isil ) so that in the rotating wave approximation the transformed 
Hamiltonian becomes time-independent. For the description of the semiconductor dynamics one can use the SBE 
accounting for degenerate valence bands J^i^Sii^ Having the dynamical equations at hand one can work out the short 
time approximation similarly to the one shown in the Introduction. However, in the short time approximation, it is 
more convenient to study directly the exciton polarization. In the coherent regime it is defined as = 
where \tpit)) is the state of the semiconductor, fi denotes the whole set of relevant quantum numbers describing the 
specific exciton (bound or unbound) state and is the exciton annihilation operator defined as (0|i?^ = In 
terms of the electron and hole operators the exciton operator is presented as 



B^i^ j diidy 0* (x, y)cs^ (x)?;^,, (y), (5) 

where (/)p(x, y) is the exciton wave functions, Va^, (y) destroys a hole in the valence band with the spin state o-^^ at the 
point with the coordinate y and c^^ (x) destroys an electron in the conduction band with the spin state at point x. 

Neglecting the Coulomb interaction and the effect of slow internal dynamics determined by small detuning lu^ 
means the semiconductor Hamiltonian is approximated by the Hamiltonian of light-matter interaction H^. In the 
rotating wave and dipole approximations the interaction Hamiltonian has the form^ i*^^'^^ 

i/e-Et^M^M^+^KW], (6) 

where e^ {t) is a piece- wise constant function of time, it takes the values when the field is off and tfj. — /dxd^ • 
E(x)(/)* (x, x) with E(x) being the envelope of the external field when the excitation is on. Here we have introduced 
the respective matrix elements of the dipole moment, which are assumed to be independent on both the electromagnetic 
field momentum and its angle of incidencci^ Thus, in this approximation one has \ip{t)) = exp(—iHf,t)\0) (throughout 
the paper we write formulas using units with h = I) and accordingly 

^^(^) = E ^ ([^- ■ • • • • •]]) . (7) 

n 

where (...) = (0| .. . |0). We calculate this series using the observation that commuting [i?^, Bl] with we arrive at 
an operator which can be presented in a form similar to Eq. ([5]). We employ this observation introducing operators 

M(n) ^Y^J d^dy [c/("i(x,y)«t (x)ct(y) + i?(")(x, y)c,(y)i;.(x)] , (8) 



which are related to each other through repetitive commutation with H^, namely, M("+^) — [iJe, [-ffe, '"']]• With 
the "initial condition" M*^"^ = B^ the operators M^") represent the even terms in series ([7]). These terms themselves 
do not contribute into since (A/^")) — 0, which in turn follows from c|0) = and (0|c^ = 0. After commuting with 
_ff, however, they produce c-number terms A'") = ^ [TJ, M*^"^]^ resulting in 

n=0 ^ ' 

In order to simplify the expressions we assume below that the quantum well can be approximated by a 2D plane. In 
this case the exciton states are characterized by the spin states of the electron and hole constituting the exciton, the 
center of mass momentum in the plane of the well, K, and other quantum numbers, n^, so that = \a^, s^, K^, n^). 
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Using the relation between M*^""*"^^ and Af*^"^ we find 

cr' ,s' 

+C/2(x,y)£.,,,(y)£:, ,,(y) - 2^^?),, (y, x)£.,,, (x)£.,,,(y) 
i?(r'/i)(x,y) = ^ [i?(?;^(x,y)£:;,,(x)f.v,,(x) 

(x,y)f:,,,(y)f.,,,,(y) - 2t/^?Uy, x)^.' (x)^:, ,(y) 

where £(j,s(x) = dcsEg^ai^)- In terms of the kernels of the operators M*^"^ the coefficients A^"^ are expressed as 
= / dx (x) with 



(10) 



' (x) - E [^;.(x)C/tnx,x) - £.,,(x)7^W(x,x)] . (11) 



Some general results can be obtained directly from Eqs. pU|) and pT|) . First, the polarization of dark excitons 
(states with helicity and ±2) is zero. Second, if the excitation pulse is circularly polarized then summation over the 
spin indices in Eqs. (|10p and (jlip reduces to the single terms with the electron and hole spins determined by and 
cr^, respectively. If the excitation pulse is linearly polarized then only the electron spin is fixed to and as a result 
different electron-hole states get coupledi^ The summation over the hole spins, as will be explicitly demonstrated 
below, to be limited to such values that meet the condition ct — (t^ G {2,0,-2}. For example, if cr^ = 3/2 (heavy 
holes) then cr can take values 3/2 and —1/2 (light hole). 

It follows from Eqs. ^ and ^ that A("+i)(x) = (x)A(") (x), where 

\£-.sM)\' (12) 

with CT^ being the second value satisfying the spin selection rule discussed above. Applying these results to Eq. ^ 
we find the general representation, which is valid for an excitation pulse with an arbitrary spatial profile in the plane 
of the quantum well, 

P^{t) =-iJ dx^^£:^(x)0;(x,x)sin[f7^(x)i] . (13) 

This result can be generalized for the more general case when the envelope amplitude of the external field is not a 
piece-wise constant and can be presented as E(x, i) = /(t)E(x). We note that the Hamiltonians H'^{t) = f{t)He 
taken at different instants ti and ^2 commute with each other. Therefore, the line of arguments used above can be 
repeated with the substitution 

t 

t^fdt'fit'). (14) 







It can be easily seen that the expression for the exciton polarization Pfj,{t) remains essentially the same 

-P^W = J '^x ^ ^^ g^(x,t)./)*(x,x)sin[n7^(x,t)] , (15) 

where f^(x, t) and fi^(x, t) are defined by the same expressions but with the time-dependent envelope function of the 
external field, i.e. they differ by the factor f{t). The phases 

t 

w^{x,t) = J dt'n^{x,t) (16) 



have the meaning of the pulse areas. Since such generalization does not bring new physics but complicates the 
discussion of the time dependence in what follows we will consider only the case of piece-wise constant amplitude of 
the external excitation. 
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FIG. 1: The ratio of two Rabi frequencies ^h+/^h- as a function of ellipticity of the excitation pulse E-^/E-. The ratio varies 
from \di/dh\ to \dh/di\. For calculations the ratio of heavy- and light-hole dipole moments is taken du/di — \/3- 



A. Single-pulse response 

First we briefly discuss the case when the excitation is a single plane wave E(x) — [Ej^e^ 4- i?_e_)e'^° ''. In this 
case 0^(x) does not depend on the coordinate and the integration in Eq. (fT3|) can be easily performed taking into 
account that 0* (x,y) = e~^^^'^^0*(x — y), where is the coordinate of the center of mass of the exciton and 
(^* (x — y) is its relative wave function. For example, for {h-\-, Is) exciton polarization we obtain 

n+,K(t) = -i^(K - Ko)d„+S+0,:+(O)-^ sminn+t), (17) 

where 0j!j_|^(O) is the value of the exciton wave function at the origin, which in the 2D approximation is 0^^(O) = 
r^^-y/2/7r with rh+ being the respective exciton Bohr radius, and 

= 2^\dH+E+\' + \di^E^\^. (18) 

Equation (fT7|) presents the Rabi oscillations of the exciton polarization. It is interesting to note that this result 
explicitly shows the statistical origin of the Rabi oscillations. Following the same line of arguments one can show 
that the oscillations would absent and the exciton polarization would increase monotonously with time if the exciton 
operators had obeyed bosonic commutation relations or if the electrons and holes were bosons. Indeed, let, for 
example, the exciton operators be bosonic then all terms in Eq. ([7]) with n > \ turn to leaving cx t. 

Because of the coupling between different exciton states the evolution of the exciton polarization is characterized 

by two Rabi frequencies, Vlh+ ~ ^i- and Q,h- = = 2\J \dh-E-\^ + \di+E+^ . The ratio between these frequencies 

is shown in Fig. [1] as a function of E+ / E- , where for calculations we have taken'^^ dh/di = 

The effect of sharing the Rabi frequency leads to the same time dependence of the heavy- and light-hole contri- 
butions into the signal with fixed helicity if the excitation is linearly polarized. If, however, the external field has 
elliptic polarization these contributions oscillate with different frequencies leading to beatings and to a nontrivial time 
dependence of the polarization state of the signal. 



B. The operator equations of motion 



Before we proceed we would like to discuss in details the spin selection rules noticed above from the perspective of 
the operator equations of motion. For this we consider an operator of the form 

P(i;x,y) Wf^aCy) +5si,s2(i)4i Wcs2(y)+ 

7rW(t)4(x)ct(y) + 7r(2)(t)c,(y)^;.(x) + A{t)5{^ - y), 
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where the sum is taken over all spin indices. The amplitudes in this expression are chosen in such a way that F{t; x, y) 
is the Heisenberg representation of the operator F{0; x, y) defined by the light-matter Hamiltonian 

i?(i;x,y)=e*^=*i?(0;x,y)e-*^=*. (20) 



The last term in Eq. (fT9|) determines the average value ^F(i;x, y)y and for the case of the exciton polarization was 
considered above. For the purpose of our discussion it suffices to consider the dynamical equations for the time- 
dependent amplitudes gai,<72{t) and so on for the case of x = y and spatially independent external field £cr,s- From 

the Heisenberg equation of motion —iF = [H^/F] we find (see Appendix [X\i 

f , ^(2) _ ^ (1) 

n -if ,tt'-^'> -iJ^^ F* 

^o-.s ~ ^^a,s' f s' ,s + i^a' ,s9(y' ,(y^ 

where the summation over the dashed spin variables is implied. With the respective initial conditions these equations 
define the Heisenberg representation of the operator F(0;x, y). For example, if all amplitudes but tt^^I are initially 
zero, then this would be the representation of the interband polarization creation operator and so on. It should be 
noted that no restrictions on the amplitudes being "bright" or "dark" are imposed in Eq. pip . Thus this is a system 
with respect to 36 unknowns and, counting degeneracy, has as many characteristic frequencies. The actual number 
of different frequencies, however, turns out to be much smaller. 

Taking the derivative of the last pair of equations with respect to time we obtain the closed system 

n^^'>--£ ,£*, ,7r'^^ - tt*^^ £*, ,£ , + 2£ -tt*^^ £, 

'^cr,s ^ '^a,s"^(j' ,s' '^a' ,s a.s''^ (J' ,s''^c ,s ^ ^<^a ,s' „i ^gi'^a' ,s, 

;-(2)__p* (2) __(2)p* ,2£* ,7r*^^ £* 

where as as well the summation over dashed spin indices is performed. Taking the complex conjugation of these 
equations we obtain the same system but with respect to (tt'^^' ,7r^^) ). Because of this symmetry the solutions of 
Eq. (|22p are divided into two classes tt^^^^ = tt*^^^ and , which correspond to A and B irreducible 

representations of the group Z2, respectively. Additional simplifications come from the structure of the coupling 
induced by the external field. Having in mind successive application for the dynamics determined by the semiconductor 
Bloch equation we consider the additional symmetries on a general ground. Due to the electric dipole selection rules 
only four elements of £„ s are non-zero. 



■£h- 








£1- 


£1+ 





, 


£h+ 



V £h+) 

where the columns are enumerated by the electron spin projections 1/2, —1/2 and rows denote different hole spins 
3/2, . . . , —3/2. We denote by Vg the vector space containing the elements of form but without any relation 
between the nonzero elements, that is the 4-dimensional subspace of the interband polarizations. For two elements m, 
w G Vf we consider / = v^u, where the composition is the usual matrix multiplication of two (rectangular) matrices 

fsi.S2 ,s^U„: ^S2- (24) 

a' 

The space spanned by such elements will be called Ve- This is a 2-dimensional subspace of the 4-dimensional space 
of the electron densities. Using Eq. ((23)) in p4|) one can see that in the basis of the electron-hole spins Ve consists 
of diagonal 2x2 matrices. Finally, for the elements u and v we consider g = vvF , which are the matrices with the 
matrix elements 



9cri,<J2 



^Va^s'Ua2,s'- (25) 



We will call the respective linear vector space V/j, which, as can be easily checked, is 8-dimensional subspace of the 
16-dimensional space of the hole-hole correlation functions. 
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The elements of Ve and Vh can be considered as operators acting on the electron and the hole spin variables, 
respectively. The important properties of these operators is that their action is reducible, namely, Vg leaves the spin 
state intact and Vh mixes only states within the groups {3/2,-1/2} and {—3/2,1/2} but not from the different 
groups. 

Using these results we obtain the modes of system p2|). These modes arc naturally classified by the 
invariant spin classes, which do not mix during the evolution, and the irreducible representations of the 
group Z2. Inspecting Eqs. one can see that the spin classes in the spin space of the electron-hole 

polarizations are Ci = {(-3/2, -1/2), (1/2, -1/2)} or {h+,l-}, C2 = {(3/2, 1/2), (-1/2, 1/2)} and C3 = 
{(3/2, -1/2), (1/2, 1/2), (-1/2, -1/2), (-3/2, 1/2)}. The first two classes correspond to the bright exciton operators, 
the last class contains the dark excitons. 



1. The bright exciton classes Ci,2 

In each class both states have the same projection of the electron spin, s = 1/2 for Ci and s — —1/2 for C2, while 
the hole spin can take two values. If we need to distinguish these values we will denote them u and u. In order to find 
the modes one can present 7ri^l(t) = tt^^s sin(wt) in Eq. (|22|). Studying the resulting system of homogeneous equations 
one finds the frequencies and (non-normalized) eigen-vectors of the modes 

where ^1 = |fcr,sP + |^5-,sP ( 4/3^ coincides with the frequency defined in Eq. and A and B denote the even and 

odd irreducible representations, respectively. Using these solutions in the first two equations of Eqs. ((2T|) we find the 
electron and hole densities coupled to the respective interband polarization operators. One can immediately see that 
these densities belong to Vg and V;i, respectively. 

As immediately follows from Eqs. pT|) and (|A2[) only _B-solutions can contribute into the average value of the 
exciton polarization, since for the symmetric solutions both terms under the integral in Eq. (|A2p cancel each other. 

Moreover, for the i3-solutions with the frequencies uP' — 01 the sum 'n'^a}s£a,s + t^'s^s^s.s vanishes, which leaves only 
solution with llP' = 4/3^ in accordance with Eq. (fT^l) . 



1 



iA,B) 



(26) 



2. The dark exciton class C3 

The dark exciton operators are specified by the amplitudes in the complement to Vs. The states are naturally 
enumerated by the hole spin state a and for each amplitude tt^.s there is a correspondent component of the excitation 
field obtained by inversion of the electron spin £a-,s with s = — s. The frequencies and non-zero components of different 
modes are found to be 

'^'^—1^1/2^ '"'3/2.-1/2^^7 ' ""-1/2,-1/2^7 ' (^'^) 

' ■ ' ^3/24/2 ' ' '^-1/2,1/2 

W^=/3!-l/2, '^l/2.1/2=~7 ' ■^-3/2,1/2 = 7 ' iAB) 

' ' ' ^^1/2, -1/2 '^•-3/2, -1/2 (^1) 

^ (/3l/2±/?_l/2)^ 7r^y2,-l/2 = ^^-^^2,-1/2 = T/3l/2, 'ri^]/2,l/2 = '^1/2,1/2 = T/3-1/2, {A), 



'3/2,-1/2 ^ "-1/2,-1/2 rHL/2i "-3/2,1/2 ^ "1/2,1/2 
= {Pl/2 ± ^-1/2)^ , 4/2,-1/2 = ^^-1/2,-1/2 = ±/3l/2, 71'-i/2,l/2 = ^^1/2,1/2 = ±/3-l/2- (B) 



Since for each tTq-.s in this class the respective components of the external field are zero, i.e. £a;s ~ 0, the dark exciton 
amplitudes, according to Eq. (|A2[) . do not contribute to A{t). Thus the average value of the dark exciton polarization 
equals to its initial value Fcr,s(i) = Fa.s{0), where Fcr,s(0) = wjcj with a and s from C3. If the initial state of the 
system is vacuum then (i^o- ^(i)) = confirming the result about absent polarization of dark excitons. 
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3. The operator dynamics 

The frequencies found = 4/3^, /J^, (/3s ± Ps)^ define all the frequencies introduced into the dynamics due to the 
interaction with the external field as defined by Eq. (|2ip . The remaining 8 frequencies can be shown to be 0. Four 
solutions (two for each class Ci and C2) corresponding to the zero frequency are presented in Eqs. (PSI) and the 
remaining 4 can be found from the first two equations of Eqs. ((2T|) . Since Eqs. (|2T|) have the form of a Schrodinger 
equation with a hermitian Hamiltonian the solutions corresponding to a; = are not secular solutions but rather 
invariants, that is they correspond to operators, which are invariant with respect to action of the external field. For 
example, as follows from Eqs. (P^ . for the exciton polarization such operators are J2a s ^'^^'^''^a^l + ^■'^■^ where the 
spin summation is restricted to a particular class Ci or C2. 

In the short-time limit only frequencies 4/3^ contribute to the Rabi oscillations of the exciton polarization. However, 
as will be demonstrated later the slow dynamics leads also to admixture of other frequencies. 



C. Multi-wave mixing response 



A more complex situation arises when the excitation pulse consists of two plane waves with non-coUinear wave- 
vectors, E(K) = E(^^(5(K — Ki) + E(^M(K — K2). This corresponds to the standard problem of excitation by two 
pulses with zero delay between them. Factoring out the exponential factor corresponding to the first pulse we can 
rewrite Eq. as 



(28) 



where F{x) is a periodic function F(x-|- AK/(27rAi^^)) = F{x.) with AK = K2 — Ki. Thus, integration over x yields 
the exciton polarization in the form 



E ^i"ni)'5(K,, - Ki - mAK) 



(29) 



where Pll™\t) are the amplitudes of multi-wave mixing (MWM) polarizations. The amplitudes Plf\t) and Pjl^^t) 
correspond to the polarization along the directions of linear response, Ki and K2, respectively. For m = —1,2 one has 
four-wave mixing polarization, m = —2, 3 correspond to six-wave mixing and so on. Performing the Fourier transform 
of the periodic factor F(x) we derive the integral representation for the amplitudes of MWM polarization 



2tt 



2tt 



(1) 



sm I 



(30) 



where = d^E^jl^ and the frequencies 57}/''^'' are defined as the amplitudes of r2p(x) corresponding to different 



,(1,2) 



pulses, =2^EJfB.P- 

In Fig. [2] we show the numerical evaluations of Eq. ([50)1 for different ratios between the amplitudes of the plane 
waves constituting the excitation along one of the forward and four-wave mixing directions. The initial stage of the 
dynamics for t — s- can be obtained directly from series © in the form (for m > 0) 



(2m -1)! r^^ ^'^'^ ^ 



(31) 



where rj^ - 

The evolution of the exciton polarization in the opposite limit t — s- cx) is found from integral representation (I30p 
using the stationary phase method 



plrHt) 



sign(4i) 



|c(l) c(2) 



: sm 
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(a) (b) 

FIG. 2: (Color online) The exciton polarization in the case when the excitation consists of two plane waves with non-coUinear 
wave vectors. The solid and dashed lines show the response in the forward (m = 1) and the four-wave mixing (m = 2) directions 
for (a) ^{.^Vi^M ' = 1> (b) = 2. The effective time is defined as t = Q^i^h. 



It shows that in this Umit the P^i™'' (t) demonstrates power decay cx t and is essentially a superposition of two 



harmonics with frequencies 



dp ± n 



(2) 



Of course, the limit t oo should not be taken too literally, since in this 



limit the main assumption w^ti ^ 1 is no longer valid and the effect of desynchronization of the MWM polarization 



of high orders (see below) becomes important. However, if fl^}'' and are not too close then after just a few 
oscillations the dependence Pli^^\t) is satisfactory approximated by its asymptotic form. For example, the complex 
behavior of the exciton polarization shown in Fig. ^Bp is a result of existence of two harmonics with the frequencies 



l(2) 



(2) 



, which exist if 51 



(1) 



If n 



(1) _ o(2) 



ri)t = VLf^ then the slow harmonic vanishes. 



Indeed, in this case, as follows from Eq. (|5D|). MWM 



polarizations can be presented in the form 



P(")(i) 



*(-i)"v;(o) 



c.(2) 



r(2) 



dt' J2m{2VLf,t') 



(33) 



where Jm{t) are the Bessel functions of the first kind. This solution asymptotically oscillates with the frequency 251^, 
which is consistent with the simple picture of the Rabi frequency renormalized due to doubling the external field. 

The interesting difference between the Rabi oscillations in the case of a single pulse excitation and the oscillations 
of the MWM polarization is that the latter asymptotically decays in time oc t^^/"^ . It can be understood taking into 
account the dynamical origin of formation the MWM polarization. As follows from Eq. (|30p the MWM polarizations 
P^™'' (t) satisfy the equation of motion of a classical tight-binding model 



(t)^f7^P^™)(<)+r^2p(™+i)(t) 



. „2 p(m-l] 



(34) 



where fi^ = 

j(m) 



dff' . Equations are supplemented by the initial conditions Pll"^'{0) = for all m and 

P^"-'(0) = for all TO except to = 0, 1, for which one has pIi°\o) = -if (^^^(0) and P^^\o) = -iS'^^UliO). This 
representation allows one to estimate the effect of the dependence of the exciton energy on the in-plane momentum. 
It leads to modification of the "eigenfrequency" , il^, by the term ~ [mAKY /2m^. Even if {AK)"^ /2m^ ^ f2^ this 
contribution becomes essential for MWM exciton polarizations of significantly high orders. However, for the most 
important case of MWM polarizations of low orders on the time scale of few Rabi flops the effect of the correction 
due to the exciton dispersion can be neglected. 

It follows from the representation (|34p that initially the excitation is localized at the sites to = 0, 1 (the directions 
of linear response). With time the energy spreads along the chain thanks to the coupling oc rf between different sites. 
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In terms of the MWM polarization the dynamics of this spreading corresponds to excitation MWM polarizations of 
higher and higher orders. The dynamical model described by Eq. ([M)) has the first integral 



(35) 



From this perspective it is clear that the decay of the Rabi oscillations is the consequence of spreading the excitation 
among MWM polarizations of different orders. 



III. THE EFFECT OF INTERNAL DYNAMICS AND THE COULOMB INTERACTION 



In the previous section we have studied in details the dynamics of the exciton polarization in the short time limit, 
when the main contribution into dynamics comes from the interaction with the external field. Here we consider 
the effect of the internal semiconductor dynamics and study how the results of the previous section appear in the 
framework of a more general description. We approach this problem describing the dynamics of the interband 
polarizations Po-,s(xi,X2) — (cs(x2)vct(xi)), the electron-electron esi,s2 (xi, X2) = (cj^ (xi)cs2 (X2)) and the hole-hole 
/icrj^o-2(xi,X2) = (wj.^ (xi)uo-2 (X2)) correlation functions by the semiconductor Bloch equation (see Appendix [B|) . 

First we notice that the most significant contribution of the full dynamics would be the production of the response, 
which was absent in the approximation used in the previous Section. This would be generating the interband polar- 
ization corresponding to dark states (class C3). It can be proven, however, that the dynamics described by Eqs. (jB5[) 
does not support excitation of the dark states. Indeed, integrating both sides of the dynamical equations with respect 
to time in the interval (0, At) with At ^ one can see that initially the spin state of the interband polarization 
belongs to Vs- Next, considering the integral over the interval {t,t + At) one can see that if at instant t the state of 
the system, i.e. p, h and e, belongs to Vg, Vh and Ve, respectively, then so does the state at t -I- At. It should be 
noted, thereby, that in quantum wells the valence band mixing does not lead to a violation of invariance of the spin 
classes since the Hamiltonian describing the mixing is an element of Vh- Thus, if initially the system is in the ground 
state then the interband polarization corresponding only to the bright states will be produced at least as long as the 
SBE and the approximation of angular independence of the dipole moment hold. 

As the next step we consider the limit of very short time response such that following the analysis of simple model ([3|) 
the contribution of the kinetic energy terms can be neglected, while other terms including non-linear are kept. The 
solution found in the previous section, Eq. (jl3p . suggests the interband polarization in the form p(x, y) cx (5(x — y). 
Substituting this ansatz into Eqs. (|B5|) with neglected kinetic energy terms and introducing a cut-off of the Coulomb 
potential such that V{0) = Vb < 00 one arrives at the system of equations independent of Vb, which is structurally 
similar to Eqs. (|A2[) and which yields the same result for the exciton polarization as Eq. p^ . The cancelation of the 
terms proportional to the Coulomb potential is particularly evident from the representation of the SBE in form (|B3[) . 

This consideration shows that the most significant correction to the approximate solution results from the kinetic 
energy terms rather than from the Coulomb interaction as one might expect considering that the Coulomb interaction 
is responsible for the non-linear terms in the SBE. 

We would like to emphasize at this point that the quantity of our main interest is the exciton polarization P^^ since 
these are the excitons that constitute true single-particle semiconductor states. As follows Eq. ([5]), is related to 
the interband polarization pa-,s through 

= y"dxdy0*(x,y)p^^,^^(y,x). (36) 

For bound exciton states this naturally introduces a characteristic spatial scale, the exciton Bohr radius rs- This 
circumstance together with the assumption that the typical spatial variation of the external excitation in the plane 
of the quantum well is small comparing to the Bohr radius Ki_2rB <C 1 allows one to simplify the SBE using the 
parametric approximation. For this we introduce new spatial variable R = (x + y)/2 and r = x — y. Next we notice 
that in the case of single pulse or two pulse excitations the solutions of the SBE are invariant with respect to either 
infinitesimal or finite (by AK/(27rAiir^)) translations in R plane. Having in mind the consecutive convolution of the 
solutions with the exciton wave function we can neglect the terms oc AKrs and leave only the dependence on R as 
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a parameter i^ii^i^ Thus, assuming for simphcity that Ki = we approximate Eqs. (|B5[) by 

ip<,,,(r,R) =i^,,,b]+£:,,,(R)(5(r)- / dr' [£,,,,(r', R)/i,,,,(r' - r, R) - R)e.',.(r' + r, R) 



-ies„s,ir,-R)^Ks,^sAe]+ J dr' [f^,^,^ (r', R)p,,,,, (r' + r, R) (r', R)p;,^,^ (r' - r, R)] , 
-z/i,„,,(r,R) =K,,^,,[h] + f dr' [5;^^,,(r', RK,,,,(-r - r', R) (r', R^^^,, (r + r', R) 



(37) 



Here £*^(r,R) = £o..s(R-)i5(r) — V{^)Pa.s{r,'R.) and the integro-difFerential operators in Eqs. jBTj are substituted by 



Hai ,cr' fa 



1 

2mcri,cr' dr"' 



fs, (r, R) - y dr' l/(r')e., (r', R)U (r - r', R), 
U'ir,-R)+H,,U,{r,-R)- f dr' V{r')K,,,,{r' ,-R)U,{r - r' ,-R). 



(38) 



System ([57]) in the short time limit when all terms in the right-hand side but those depending on the external field are 
neglected reproduces the results obtained in the previous section from Eq. ([7]). In particular the multi-wave mixing 
response is obtained by the Fourier transform over the parameter R similarly to Eq. pop . In what follows we will 
use Eqs. p7|) in even simpler form assuming the single-pulse excitation only and thus omitting the dependence on R. 
It suffices for our purpose of discussion of the effect of the internal dynamics and the Coulomb interaction. We also 
omit the off-diagonal elements of the holes Hamiltonian H^^^a^ describing the mixing of the valence bands. Thus we 
characterize these bands by ruhh and mih that is by the masses of the heavy- and light-holes, respectively. 

Considering the terms besides the kinetic energy in the right-hand side of these equations as sources the responses are 
governed by the free-particle propagators for the correlation functions and the Coulomb propagator for the interband 
polarization. The fundamental property of these propagators is that in the short time limit i — s- they turn to spatial 
(5-functions owing to fast spatial oscillations of the kernel 



if(ri,r2;t) 



— exp [jTOlri — r2p/2il 

, J- L J 



(39) 



In order to evaluate the parameter m we take into account that the evolution of the semiconductor initially being 
in the ground state is spanned by the spin classes Ci^2- Thus for the electron correlation function m = oo so that 
the propagator remains local. For the hole correlation functions this parameter is finite for the off-diagonal elements 



(to 



-1 

Ih 



^hh 



) . For the interband polarization in the short time limit the Coulomb propagator can be 



approximated by the free-particle propagator-^^*^ with to « rrixh = {nT'e^ + 'm'hh)~^ and m w ruxi = (toJ^ -|- mihY 
where rrixa are the excitons' reduced masses and TOe is the electron mass. The effect of the dispersion of the electrons 
and holes on the exciton polarization is expressed as a distortion of the exciton wave function under the action of 
the free-particle propagators. The time scale, at which the distortion becomes essential, can be estimated as the time 
required for the initial (5-shape of the kernel to acquire the width of the order of the exciton Bohr radius. This leads 
to the estimate 



47r 



(40) 



The numerical value of the typical time scale for GaAs is determined by mhh = 0.45too, mih = 0.082too, m = toa 
and tb = r^h = / 'i-'nixh^^ , where mo is the electron mass in empty space and ef, = 13 is the background dielectric 
function. Using these values in Eq. (|40|) we find tc ~ 20 fs. Thus for pulses with duration shorter than 20 fs and the 
intensity sufficiently high to produce several Rabi flops, the approximation used in the previous Section is satisfactory. 
However, as will be shown below the exciton polarizations corresponding to the states in the discrete spectrum are 
less sensitive to the deviation of the dynamics from its short time limit. 

First we discuss the general effect of the electron-hole dispersion on the Rabi oscillations. For this we consider 
the case of excitation by a single circularly polarized pulse of high intensity. For such excitation we neglect the 
contribution of the non-linear Coulomb terms (see below) and additionally simplify the internal dynamics assuming 
the effects of the detuning and the energy offsets to be small at the time scales under consideration. Since in this 
case the translational invariance in the plane of the quantum well is restored it is more convenient to solve Eqs. ([37]) 
rewriting them in the momentum representationi^iiSiiH The Laplace transformation p(A) = p{t)e~^*dt of the 
solution is found to be 
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FIG. 3: The effect of the electron-hole dispersion on the Rabi oscillations. The time dependence of the real part of the exciton 
polarization is shown for flcr,s/i^iT.3{l/rB) = 5 (solid line) and QcT.s/i-Ocr,s{l/rB) ~ 10 (dashed line). 



where a;o-,s(k) = k^/2mxa and fl^ g — 4(|£ct.sP + Ii^-ct.sP) with a being the hole spin state comphmentary to a in the 
respective spin class. In the case under consideration because of circular polarization of the excitation pulse only 
interband polarization corresponding to the same helicity is created, thus, if £^,s ^ then E^^s ~ and vice versa. 
We, however, write down the "full" expression for ilg-.s in order to show the relation with Eq. (I12p . The poles of the 
right-hand side of Eq. (1411) determine the characteristic frequencies yielding the time dependence in the form 

p.^.(t,k) = -^ f^'- jsin[f7.^,(k)t] - 2z^^sin^[17.,,(k)t/2]) 

e,,.(t,k) = KAt.^) ^ ,^ ti?^',. , sin^[r2.,.(k)t/2] 

where r2o-,s(k) — ^J^l^s + ^o-,s(k)^. It is interesting to note the similarity to the simple dynamical models discussed 
in Introduction [see Eq. ([4])]. 

The exciton polarization is found by convoluting Eq. (j42p with the exciton wave function 0o^^s(k) = 2^/2TTr„^s{l + 
fc^r^^)~^/^ according to Eq. In Fig. [3] we plot the time dependence of the exciton polarization for different 

intensities of the external excitation. Since the exciton wave function mostly localized inside the region k ~ l/r^ and 
by assumption ujn^s(l/rB) < i^o-.s this convolution does not lead to appearance of new frequencies, thus the oscillations 
of the exciton polarization are characterized by a slightly modified Rabi frequency Oo-,s- 

There are two significant corrections to the picture of the oscillations obtained in the previous section. First, there 
is slower than exponential decay due to dephasing of the components corresponding to different k. The decay can be 
estimated considering the long time asymptotic t — > cx). Since the phase in Eq. (|42p reaches its stationary points at 
the ends of the k interval, i.e. at fc = 0, it implies the asymptotic power decay of the amplitude as 1/t independent 
on the intensity of the excitation. The second modification is related to the second term in the brackets in Eq. ([1^ 
[Fig. Od]. This term has the form of oscillations with the Rabi frequency Qa,s near a constant level, which position 
strongly depends on the intensity of the excitation for fla-,s ^ '^^it.s{^/'''b) and slowly decays towards for higher 
values of the Rabi frequency as shown in Fig. [H 

Solutions allow us to estimate the effect of the Coulomb interaction and to show that the contribution of the 
non-linear terms into the right-hand side of equations of motion (|37p vanishes with increasing the intensity of the 
excitation. Physically this is supported by, first, observation that the interband polarization and the single particle 
correlation functions as functions of the external excitation are limited from above and, hence, the nonlinear terms 
are scaled as the exciton binding energy^^ , while the other terms increase with increasing the intensity of the external 
field. Second, one can also apply the result obtained from Eq. (|B3I) about vanishing the non-linear terms in the 
limit of (5- functional profile of the interband polarization and the single-particle correlation functions. As follows from 
Eqs. (|42p with increasing the intensity of the excitation the spatial distribution of Pa.s, ha-^a- and e^.s narrows with the 
characteristic width oc 1/kfi, where fcn — ^y2/rn^^^^^l^ . 
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FIG. 4: The dependence of the constant level on the Rabi frequency Q,a,s/ijJtT,s{'^/rB)- 



There are different ways to estimate the contribution of the non-linear Coulomb terms into the dynamics. In order 
to illustrate the application of the physical ideas mentioned above we consider the contribution into the right-hand 
side of the equation of motion for the interband polarization coming from the term written in the second line in 
Eq. (jB3[) . In the simplest case of a single pulse excitation this term yields an addition to the equation with respect 
to Pcr.s(i,k) in the form 

A5,,,(k) = j dk'F(k - k') [p,,,(k')n(k) - p,,,(k)n(k')] , (43) 

where ?t.(x) = es,s{'^)-\-ha^cr{'^) and we have taken into account that the dynamics is restricted to a particular spin class. 
We illustrate the typical line of calculations leaving only the first term in the expression for po-_s(k) and approximating 
the oscillating terms by 1 thus presenting 

A5,,,(k) ^ -1^$^ [ dk' ^J^^^Z^'^^ ^ In^A^) - r2.,.(k')l . (44) 



Rescaling the wave vectors and using for the single-pulse excitation fla^s ~ S^^.s we obtain 

1 / ET-rr . / k 



AS,A^)^-^-—JJ-^C^ (45) 



where exa = '^rUxaS^/^l is the exciton binding energy and 



^^^^ - 1 JVT^WT^JW^y ^''^ 

The dependence of this function on k is shown in Fig.[5l It varies slowly over the interval k < 1/2, where its magnitude 
is close to C(0) = 7r^/V2 — 8r^(5/4)^/7r « —4.67. Thus, if the Rabi frequency exceeds the exciton binding energy or, 
equivalently, if the typical spatial scale of the interband polarization l/Zcji is smaller than the exciton Bohr radius 
the effect of the non-linear terms on the time dependence of the exciton polarization reduces to a relatively small 
modification of the amplitude of the external field. 

This argument allows one to reformulate the criterion of validity of the short time approximation for the exciton 
polarization in terms of relation between the exciton Bohr radius tb and the typical spatial width of the interband 
polarization due to the interaction with the external field ~ 1/kn. If rskfi 3> 1 then the result of convolution of 
the exciton wave function with the interband polarization pa-^s does not depend on the details of the dependence of 
Pcr,s(k) on k. Indeed, since in this case the main contribution into the integral P ^ J rfx?/;*(x)p(x) comes from the 
small area near x = we can approximate the integral by 

P « V*(0) / dxp(x) = (27r)2^/>*(0)p(k = 0) (47) 



which turns Eq. (j42p into Eq. (|17p exactly (including the cancelation of the factor (27r)^). This consideration shows 
that the exciton polarization corresponding to the states in the discrete spectrum is less sensitive to the deviations 
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FIG. 5: The dependence on k = fc/fcn of the integral in Eq. (146 |l describing the effective modification of the source due to the 
nonlinear terms in the limit of high intensity of the external field. 



from the short-time dynamics of the system than the interband polarization. The latter contains information about 
all exciton states including those from the continuum spectrum whose spatial variation is no longer characterized by 
the exciton Bohr radius. The condition rsfcn = ^/IP'^.s/ ^xa) 3> 1 imposes the limitation on the intensity of the 
external field that the period of the Rabi oscillations should be smaller than ~ 200 fs (for GaAs). 

An interesting effect of the dispersion is revealed when the excitation pulse is not circularly polarized. In this case 
both states constituting a spin class are involved in the dynamics resulting in the appearance of a new characteristic 
frequency. In order to illustrate this effect we consider the case when the intensity of the excitation is significantly 
high so that the effect dispersion of the hole-hole correlation function, which is quantified by the parameter = 
,j, can be treated perturbatively. In order to simplify the expression we assume additionally that both 



''hh 



heavy- and light-hole excitons are described by the same dispersion law uja^si^)- Under these assumptions the Laplace 
transform of the interband polarization can be written as 



(A,k) 



1 



27r2A 



[A2 



(k)][(A- 



.(k)/2)2 



(k)/4] 



(48) 



where pal{\,ii) is given by Eq. (|¥T1) and iLjg-,a-(k-) — k^/'^.mg-a- with m^]j. — m^^ — i-e. Im-s-.o-l = m-A- The 

perturbation term in the right-hand side of Eq. (j48|) has poles at the same Rabi frequency il^ ^ (k) as in the case of 

the circularly polarized excitation and at the frequency ~ ^ (k) /4 yielding a component of the exciton polarization 
oscillating with doubled period. We would like to emphasize that the new frequency was initially present among 
the eigenfrequencies of the bright exciton classes [see Eq. d^^ ]. Its appearance, however, in the time evolution 



of the exciton polarization was prohibited in the dispersionless limit by the symmetry p' 



(2)* ^ of dynamical 

equations (|22p. The exciton finite masses break the symmetry allowing for the frequency il^j,(k)/4 to contribute to 
the evolution. With this regard it could be speculated that if the dark excitons would be involved in the dynamics 
their presence could be traced by admixture of the frequencies specific for the dark exciton spin class C3 [see Eq. (P7|) ]. 

It should be emphasized that neither effect of decay of the Rabi oscillations nor the admixture of "extra" frequencies 
are present in atomic analogues of the semiconductor Rabi oscillations. The decay is essentially caused by the spatial 
spread of the true (single-particle) eigenstates of the system. The appearance of additional frequencies is related to 
the existence of different kinds of holes with continuous spectra characterized by different masses. 



IV. CONCLUSION 



We have studied the short time semiconductor response with respect to resonant high intensity excitations. The 
main results are obtained neglecting the slow dynamics, that is on the time scale provided by detuning. We have 
calculated exactly the excitation polarization — (?A(t)| in the coherent limit. We have found that if the 

semiconductor is excited by a field consisting of one or two plane waves the polarization of bright excitons demonstrates 
Rabi oscillations, while the polarization of dark excitons identically zero. The main difference between the single-wave 
and two- wave excitations is that in the latter case the Rabi oscillations of the multi-wave mixing (MWM) polarizations 
are generally a superposition of two harmonics (Oi > 5^2) unless the amplitudes of the waves satisfy a special relation. 
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For example, if both waves are circularly polarized with the same helicity then J7i — f22 ~ 2dE, where E is amplitude 
of the weaker wave, and the slow harmonic ^2 collapses only if the amplitudes of the both waves are the same. Also, 
in the case of the two pulse excitation oscillations of the exciton polarization decay cx t^^/^ owing to redistribution 
of the excitation over the multi-wave mixing polarizations. 

In order to establish the dynamical origin of the vanishing polarization of the dark states we considered the operator 
dynamics governed by the external field. We found that there exist three invariant spin classes, which do not mix 
with the evolution of the system. Two of these classes correspond to bright excitons and the third one contains all 
the dark states. It has turned out that the operator dynamics is described by six frequencies. Only two of those 
frequencies show up as the Rabi frequencies in the time evolution of the exciton polarizations (one for each bright 
classes), while others are prohibited by a hidden symmetry. 

The short time approximation used in the first part of the paper appears naturally in the more general context of 
the semiconductor Bloch equation (SBE). We use the SBE in order to discuss the effect of the Coulomb interaction 
and the effect of the continuous spectra of the electrons and holes on the results obtained in the limit of immediate 
response with respect to intensive excitation. We found that if the system is initially in the ground state then the 
SBE preserves the invariant spin classes. Thus, it proves that as long as the SBE holds and the approximation of 
angular independent dipole moment is justified the dark excitons do not appear in the dynamics. 

An interesting result of studying the short time approximation within the context of the SBE is that the Coulomb 
interaction does not lead to essential changes if the intensity of excitation is sufficiently high compared to the exciton 
binding energy. Moreover, including the nonlinear exchange terms directly to the equations of the short time limit 
(i.e. without the kinetic energy) adds nothing new since they cancel each other. In turn, the effects of the continuous 
spectrum of the electrons and holes are much more significant. We show that the electron-hole dispersion while does 
not change significantly the frequency of the oscillations of the exciton polarization leads to the power law decay of 
the oscillations. An interesting effect is related to the mass difference between heavy- and light-holes. It breaks the 
symmetry that selects only two Rabi frequencies out of six that are introduced into the system by the external field 
and leads to admixture of the second nontrivial frequency for a specific bright class. 
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APPENDIX A: THE HEISENBERG REPRESENTATION INDUCED BY THE EXTERNAL FIELD 



In order to derive Eqs. (j2ip one can either consider a regularization of the ^-function, which appears in the anti- 
commutation relations and in the right-hand-side of Eq. or consider a more general form of Eq. with the 
amplitudes being scalar functions, which are convoluted with the electron and hole operators, 



F{t) = A{t) +J2 J ^^'^y -9-1 W"-^ (y) 



(Al) 



Furthermore the Heisenberg equation of motion produces the dynamical equations with respect to the amplitudes 

Ai,s2(x,y) =if<x',si(x)7r^^|^^(x,y) - j4!^^^(y,x)f;,^^^(y). 



<Ti ,<T2 (x, y) = {-^)T^^als' (y . x) - «7r^^^|,, (x, y)£l^_ (y) , 

'Tl!](x,y) = - if^,s'(x)/s,s'(y,x) - i5iff,^'(x,y)f^/,5(y), 
7T(2)(x,y) (x)/,,,,(x,y) +zg,,,,(y,x)£;, ,(y). 



(A2) 



A=i 



T^i'is' (X, X)£:;, (x) - -K^V^, (X, Ti)£„, y (X 



(2) 



where the summation over dashed spin variables is implied. The last equation gives the time dependence of (^F(t) 

(compare with Eq. pT|) ). Considering the dynamics of the diagonal, x = y, values of the amplitudes in the case of 
spatially homogeneous excitation £cr,s(x) — £a,s one obtains Eq. (PT|) . 
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It should be emphasized that Eqs. (jA2[) are not a semiconductor Bloch equations. These are equations with respect 
to amphtudes of different operators entering the Heisenberg representation of the operator F{t). The semiconductor 
Bloch equation, in turn, is written for average values of the specific operators [see e.g. Eq. (IBSp ]. While these equations 
have similar structure (in the linear approximation for the SBE) they are different as can be see comparing closely 
Eqs. (|A2|) and l|B5|) . In order to find an average value of a particular operator using Eqs. (|A2|) one needs to solve 
these equations for specific initial conditions and then use this solution in the last equation of Eqs. (|A2p . 



APPENDIX B: THE SEMICONDUCTOR BLOCH EQUATION IN THE COORDINATE 

REPRESENTATION 



In the main text we consider the semiconductor response under the action of an external field, which is spatially 
inhomogeneous in the plane of the quantum well. Because the system does not posses the translational invariance it 
is convenient to use the semiconductor Bloch equation in the coordinate representation. 

The Hamiltonian of an excited semiconductor in the rotating wave approximation has the form 



H = Hel + Hh + He + He 



(Bl) 



where Hei and H^ are the standard electron and hole single-particle Hamiltonians^ii^, He is the Hamiltonian of the 
Coulomb interaction and H^ describes the light-matter interaction 



He = j £<y,s (x)4 (x)ct (x) + h.c. 



(B2) 



where ^^.^(x) — (s|p|(t) •E(x), operator wj.(x) creates a hole with the spin state a at point x and c|(x) is the respective 
electron creating operator. The equations of motion for the interband polarizations Po-^s(xi,X2) = (cs(x2)wct(xi)), 
the electron-electron esi^s2(xi,X2) = (c]!^ (xi)cs2 (X2)) and the hole-hole /icri,CT2(xi,X2) = (t;^^ (xi)i)ct2 (X2)) are derived 
from the equation of motion d (/) / dt = i {[H, /]) using the Hartree-Fock approximation for the terms describing the 
Coulomb interaction. 

Straightforward calculations give immediately quite cumbersome system of equations. We provide explicitly only 
one equation, which shows the structure of the non-linear Coulomb terms 

«PfT,s(xi,X2) = K^^s[p] + [C^a(xi) - t/A(x2)]pCT,s(xi,X2) 

-J dx' [V^(xi-x')-V^(x2-x')][/i.',^(x',xi)p,,,,(x',X2)-p.^,'(xi,x')e,-,,(x',X2)] (B3) 

-|-£o.,s(xi)(5(xi - X2) - /lcr',o-(x2,Xi)fc,/,5(x2) - f (Xl )es' ,s (Xl , X2) . 

Here and below a summation over dashed spin indices is implied. This equation clearly demonstrates that at the 
diagonal Xi = X2 the contributions of the Coulomb terms cancel each other. In Eq. (jB3p [/a(x) describes an effective 
background potential created by the local imbalance between the electrons and holes 



f/A(x) = y"dx'V^(x-x')[e,.,,.(x',x')-/i.'.a'(x',x' 



(B4) 



For a semiconductor excited by a single plane wave this potential vanishes in an overall neutral system due to the 
translational invariance. In the two wave excitation setup the effect of Ua is small for bound exciton states if the 
order of the multi-wave mixing response is not too high. 

The dynamical equations significantly simplify after the terms being rearranged to form renormalizations of the 
single-particle energies and the coupling between the charge densities and the interband polarizations 



ipc,s(xi,X2) =K^^s[p] +fcr,s(xi)5(xi -X2) 



dx.' 



,s (x' , X2 ) /icr' ,cr (x' , Xl ) - (x' , Xl ) e^/ ,s (x' , X2 ) 



-zesi,s2(xi,X2) ^^-ft'si^ssN + J dx.' ,,^(x',Xi)p<^/,52(x',X2) - f^',s2(x',X2)p;,^^^(x',xi; 

-i/lo-i,£T2(xi,X2) + / dx' ,,,(xi,x')p^2^s/(x2,x') - ^a2,s'(x2,x')p^i,s'(x;i,x') 



(B5) 



where we have introduced the modified coupling between the interband polarization and the charge densities 

5„,,(x,y) =£,,,(x)<5(x-y)-F(x-yK,,(x,y). (B6) 
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We represent the time dependence of the external field in the form E(t) = E{t)e where the amplitude E{t) is 
assumed to be constant during the excitation. In the rotating frame the operators 



K„^s[p] =-ffcr,o-'Pcr'.s +Pcy,s'Hs'^s ^ V {yii - X2 (xi , X2 ) - f2po-,s (xi , X2 ) , 



K \e\ =H le , — e ,H, 



(B7) 



^O"! ,(7' Hfji 



are expressed in terms of the integro-differential operators H. The action of these operators is defined by 
1 32 



^(Ti.a' fa' 



7, -7^ + K(xi) + C/a(xi) + Hs, 



/si(xi,X2)- / dx' F(xi - x')esi,s'(xi,x')/s/(x',X2), 



1 



92 



■/^'(xi,X2) + [Vh(xi) - Ua{xi) + iJffJ /^i(xi,X2) 



- y"dx' t/(xi - x')/i^j^^/(xi,x')/^/(x',X2) 



(B8) 



with Ve and Vh being the electron and hole confinement potentials, which effect on the spatial motion of the electrons 
and holes is estimated by the ^-functional approximation for the quantum well. The parameters Hs and -ffo- denote the 
energy offsets of the respective bands. The diagonal elements of maimer 2 are the masses of the light- and heavy-holes, 
while the off-diagonal elements account for the valence bands mixing. 

Thus, K^^s has the meaning of a two-particle Hamiltonian with the Coulomb interaction between them, while the 
properties of free particles are described by the Hamiltonians Hf,^^„^ and -ffsi,s2- 
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